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Abstract. Havel in 1955 [28], Erdos and Gallai in 1960 [21], Hakimi in 
1962 [26], Ruskcy, Cohen, Eades and Scott in 1994 [69], Barnes and Sav- 
age in 1997 [6], Kohncrt in 2004 [49], Tripathi, Vcnugopalan and West in 
2010 [83] proposed a method to decide, whether a sequence of nonnega- 
tive integers can be the degree sequence of a simple graph. The running 
time of their algorithms is D.[n^] in worst case. In this paper we pro- 
pose a new algorithm called EGL (Erdos-Gallai Linear algorithm), whose 
worst running time is 0(n). As an application of this quick algorithm we 
computed the number of the different degree sequences of simple graphs 
for 24, . . . , 29 vertices (see [74]). 



1 Introduction 

In the practice an often appearing problem is the ranking of different objects 
as hardware or software products, cars, economical decisions, persons etc. A 

Computing Classification System 1998: G.2.2. [Graph Theory]: Subtopic - Network 
problems. 
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typical method of the ranking is the pairwise comparison of the objects, as- 
signment of points to the objects and sorting the objects according to the sums 
of the numbers of the received points. 

For example Landau [51] references to biological, Hakimi [26] to chemical, 
Kim et al. [45], Newman and Barabasi [61] to net-centric, Bozoki, Fiilop, Poesz, 
Keri, Ronyai and Temesi to economical [1, 10, 11, 42, 80], Liljeros et al. [52] 
to human applications, while Ivanyi, Khan, Lucz, Pirzada, Soter and Zhou 
[30, 31, 38, 65, 67] write on applications in sports. 

From several popular possibilities we follow the terminology and notations 
used by Paul Erdos and Tibor Gallai [21]. 

Depending from the rules of the allocation of the points there are many 
problems. In this paper we deal only with the case when the comparisons have 
two possible results: either both objects get one point, or both objects get zero 
points. In this case the results of the comparisons can be represented using 
simple graphs and the number of points gathered by the given objects are the 
degrees of the corresponding vertices. The decreasing sequence of the degrees 
is denoted by b = (bi , . . . , bn). 

From the popular problems we investigate first of all the question, how 
can we quickly decide, whether for given b does exist there a simple graph 
G whose degree sequence is b. In connection with this problem we remark 
that the main motivation for studying of this problem is the question: what is 
the complexity of deciding whether a sequence is the score sequence of some 
footbah tournament [24, 32, 35, 36, 43, 44, 54]. 

As a side effect we extended the popular data base On-line Encyclopedia of 
Integer Sequences [72] with the continuation of contained sequences. 

In connection with the similar problems we remark, that in the last years 
a lot of papers and chapters were published on the undirected graphs (for 
example [8, 9, 12, 16, 19, 29, 37, 41, 55, 68, 81, 83, 84, 85]) and also on 
directed graphs (for example [7, 11, 14, 23, 24, 30, 31, 33, 38, 45, 48, 50, 57, 
58, 63, 65, 64, 66]). 

The majority of the investigated algorithms is sequential, but there are 
parallel results too [2, 18, 20, 36, 60, 62, 77]. 

Let I, u and m integers (m > 1 and u > I). A sequence of integer numbers 
b = (bi , . . . , bra] is called (I, u, m)-bounded, if I < bi < u for i = 1 , . . . , m. 
A [l,u,m)-bounded sequence b is called [I, u,m) -regular, if bra ^ ^-m.-^ ^ 
• • • > bi. An (I, u, m)-regular sequence is called [l,u,m)-even, if the sum of 
its elements is even. A (0, n— ^ ,n) -regular sequence b is called n-graphical, if 
there exists a simple graph G whose degree sequence is b. If I = 0, u = n — 1 
and m = TL, then we use the terms n-bounded, n-regular, n-even, and n- 



232 



A. Ivanyi, L. Lucz, T. F. Mori, P. Soter 



graphical (or simply bounded, regular, even, graphical). 

In the following we deal first of all with regular sequences. In our definitions 
the bounds appear to save the testing algorithms from the checking of such 
sequences, which are obviously not graphical, therefore these bounds do not 
mean the restriction of the generality. 

The paper consists of nine parts. After the introductory Section 1 in Sec- 
tion 2 we describe the classical algorithms of the testing and reconstruction of 
degree sequences of simple graphs. Section 3 introduces several linear testing 
algorithms, then Section 4 summarizes some properties of the approximate 
algorithms. Section 5 contains the description of new precise algorithms and 
in Section 6 the running times of the classical testing algorithms are pre- 
sented. Section 7 contains enumerative results, in Section 8 we report on the 
application of the new algorithms for the computation of the number of score 
sequences of simple graphs. Finally Section 9 contains the summary of the 
results. 

Our paper [37] written in Hungarian contains further algorithms and simu- 
lation results. [35] contains a short summary on the linear Erdos-Gallai algo- 
rithm while in [36] the details of the parallel implementation of enumerating 
Erdos-Gallai algorithm are presented. 

2 Classical precise algorithms 

For a given n-regular sequence b = (bi,...,bn) the first i elements of the 
sequence we call the head of the sequence belonging to the index i, while the 
last n — i elements of the sequence we call the tail of the sequence belonging 
to the index i. 

2.1 Havel-Hakimi algorithm 

The first algorithm for the solution of the testing problem was proposed by 
Vaclav Havel Czech mathematician [28, 53]. In 1962 Louis Hakimi [26] pub- 
lished independently the same result, therefore the theorem is called today 
usually as Havel-Hakimi theorem, and the method of reconstruction is called 
Havel-Hakimi algorithm. 

Theorem 1 (Hakimi [26], Havel [28]). /f n > 3, then the n-regular sequence 
b = (bi, . . . ,bn) is n-graphical if and only if the sequence b' = [hi — 1,b3 — 
1, . . . ,bb, - l,bb,+i - 1,bbi+2, . . . ,bn) is (n- \ )-graphical. 



Proof. See [9, 26, 28, 37]. 



□ 
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If we write a recursive program based on this theorem, then according to 
the RAM model of computation its running time will be in worst case O(rL^), 
since the algorithm decreases the degrees by one, and e.g. if b = ((n — 1)^), 
then the sum of the elements of b equals to @[n^). It is worth to remark that 
the proof of the theorem is constructive, and the algorithm based on the proof 
not only tests the input in quadratic time, but also construct a corresponding 
simple graph (of course, only if it there exists). 

It is worth to remark that the algorithm was extended to directed graphs 
in which any pair of the vertices is connected with at least a > and at most 
b > a edges [30, 31]. The special case a = b = 1 was reproved in [23]. 

In 1965 Hakimi [27] gave a necessary and sufficient condition for two se- 
quences a = [a^ , . . . , an) and b = (bi , . . . , bn.) to be the in-degree sequences 
and out-degree sequence of a directed multigraph without loops. 

2.2 Erdos-Gallai algorithm 

In chronological order the next result is the necessary and sufficient theorem 
published by Paul Erdos and Tibor Gallai [21]. 

For an n-regular sequence b = (bi , . . . , bn) let = bi + • • • + bt. For given 
i the elements bi, . . . , b^ are called the head of b, belonging to i, while the 
elements b^+i , . . . , b^ are called the tail of b belonging to i. 

When we investigate the realizability of a sequence, a natural observation is 
that the degree requirement of a head is covered partially with inner and 
partially with outer degrees (with edges among the vertices of the head, resp. 
with edges, connecting a vertex of the head and a vertex of the tail). This 
observation is formalized by the following Erdos-Gallai theorem. 

Theorem 2 (Erdos, Gallai [21]) Let n > 3. The n-regular sequence b = 
(bi , . . . , bn.) is n-graphical if and only if 



n 




even 



(1) 



i=l 



and 



n 



Xbi<j(j-1)+ Y. m^O,^^) = 1,..., 



(2) 




Proof. See [9, 15, 21, 70, 83]. 



□ 
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Figure 1: Number of regular and even sequences, and the ratio of these numbers 

Although this theorem does not solve the problem of reconstruction of 
graphical sequences, the systematic application of (2) requires in worst case 
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(for example when the input sequence is graphical) 0{n ) time. 

Recently Tripathi and Vijay [83] published a constructive proof of Erdos- 
Gallai theorem and proved that their construction requires O(n^) time. 

Figure 1 shows the number of n-regular (R(rL) and n-even (E(tl) sequences 
and their ratio (E(n)/R(rL) for n = 1 , . . . , 38. According to (34) the sequence 
of these ratios tends to ^ as n tends to oo. According to Figure 1 the conver- 
gence is quick: e.g. E(20]/R(20) =0.5000006701511. 

The pseudocode of Erdos-Gallai see in [37]. 

3 Testing algorithms 

We are interested in the investigation of football sequences, where often ap- 
pears the necessity of the testing of degree sequences of simple graphs. 

A possible way to decrease the expected testing time is to use quick (linear) 
filtering algorithms which can state with a high probability, that the given 
input is not graphical, and so we need the slow precise algorithms only in the 
remaining cases. 

Now we describe a parity checking, then a binomial, and finally a headsplit- 
ting filtering algorithm. 

3.1 Parity test 

Our first test is based on the first necessary condition of Erdos-Gallai theorem. 
This test is very effective, since according to Figure 1 and Corollary 14 about 
the half of the regular sequences is odd, and our test establishes in linear time, 
that these sequences are not graphical. 

The following simple algorithm is based on (1). 

Input, n: number of the vertices (n > 1); 
b = (bi , . . . , hn): an n-regular sequence. 

Output. L: logical variable (L = False shows, that b is not graphical, while 
the meaning of the value L = True is, that the test could not decide, whether 
b is graphical or not). 

Working variable, i: cycle variable; 
H = (Hi , . . . , Hn): Hi is the sum of the first i elements of b. 

Parity-Test (n, b, L] 

01 Hi =bi 

02 for i = 2 to n 

03 Hi = Hi_i + bi 
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04 if Hn odd 

05 L = False 

06 return L 

07 L = True 

08 return L 

The running time of this algorithm is Q[n) in all cases. Figure 1 characterizes 
the efficiency of Parity-test. 

(1) is only a necessary condition, therefore Parity-Test is only an approx- 
imate (filtering) algorithm. 

3.2 Binomial test 

Our second test is based on the second necessary condition of Erdos-Gallai 
theorem. It received the given name since we estimate the number of the inner 
edges of the head of b using a binomial coefficient. Let Tt = bi^+i + • • • + b^. (i = 
l,...,n). 

Lemma 3 // n > 1 and b is an n-graphical sequence, then 

Hi<i(i-1)+Ti (i = l,...,n-l). (3) 

Proof. The left side of (3) represents the degree requirement of the head of b. 
On the right side of (3) i(l — 1 ) is an upper bound of the inner degree capacity 
of the head, while is an upper bound of the degree capacity of the tail, 
belonging to the index i. □ 

The following program is based on Lemma 3. 

Input, n: number of the vertices (n > 1); 
b = (bi , . . . , hn): an n-regular even sequence; 
H = (Hi , . . . , Hn): Hi^ the sum of the first i elements of b; 
Hq: auxiliary variable, helping to compute the elements of H. 

Output. L: logical variable (L = FALSE signals, that b is surely not graphical, 
while L = True shows, that the test could not decide, whether b is graphical). 

Working variables, i: cycle variable; 
T = (Ti , . . . , Tn): Ti the sum of the last n — i elements of b; 
Tq: auxiliary variable, helping to compute the elements of T. 

BlNOMIAL-TEST(n, b, H, L) 

01 To = 

02 for i = 1 to n - 1 
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03 Ti = Hn - Hi 

04 if Hi > t(i - 1 ) + Ti 

05 L = False 

06 return L 

07 L = True 

08 return L 

The running time of this algorithm is 0(n) in worst case, while in best case 
is only 0[^). 

According to our simulation experiments Binomial-Test is an effective 
filtering test (see Figure 2 and Figure 3). 

3.3 Splitting of the head 

We can get a better estimation of the inner capacity of the head, than the 
binomial coefficient gives in (3), if we split the head into two parts. Let [t/2J = 
hi, p the number of positive elements of b. Then the sequence (bi, . . . ,h]^^) 
is called the beginning of the head belonging to index i and the sequence 
(bht+i ) ■ ■ • 5 bi) the end of the head belonging to index i. 

Lemma 4 // n > 1 and b is an n-graphical sequence, then 

Hi < min(min(Hin, Tn — Ti, hi(n — t)) 

+ min(Hi - W^^Jn - Ti, (i - h.i)(rL - i)),Ti) 

+ min(hi(i - hi) + (^^ + '^"^ (i = 1 , . . . , u), (4) 

further 

min(Hhi,Tn-Ti,hi(n-i)) + min(Hi-Hh^,Tn-Ti, (i-h.i)(n-i)) < Ti. (5) 

Proof. Let G be a simple graph whose degree sequence is b. Then we divide 
the set of the edges of the head belonging to index i into five subsets: (Si^i ) 
contains the edges between the beginning of the head and the tail, [Si^i] the 
edges between the end of the head and the tail, Si^3 the edges between the parts 
of the head, Si^4 the edges in the beginning of the head and Si^s the edges in 
the end of the head. Let us denote the number of edges in these subsets by 
Xi,i, . . . , Xi^5. 

Xij is at most the sum H^j of the elements of the head, at most the sum 
Tn. — Ti of the elements of the tail, and at most the product h.i(n — i) of the 
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elements of the pairs formed from the tail and from the beginning of the head, 
that is 

Xi,i < min(HH^,Tn-Ti,hi(n-i)). (6) 
A similar train of thought results 

Xt,2 < min(Hi - HH,,Tn - Ti, {i - ht)(n - i)). (7) 

Xi^3 is at most h,t(i — hi) and at most H^, implying 

Xi,3 <min(hi(t-hi),Hi). (8) 

Xi^4 is at most (2^) and at most H^j, implying 

Xi,4<min(Q'^,HKj, (9) 

while Xi,_5 is at most {^ 2^) and at most Hi — Hhj, implying 

X.S <('-"'). (10) 

A requirement is also, that the tail can overrun its capacity, that is 

Xi,i + Xi,2 < Ti. (11) 

Summing of (6), (7), (8), (9), and (10) results 

Hi < Xi,i + Xi,2 + Xi,3 + 2Xi,4 + 2Xi,5. (12) 

Substituting of (6), (7), (8), (9), and (10) into (12) results (4), while (11) is 
equivalent to (5). □ 

The following algorithm executes the test based on Lemma 4. 

Input, n: the number of vertices (n > 1); 
b = (bi, . . . ,bn): an n-even sequence, accepted by Binomial- Test; 
H = (Hi , . . . , Hn): Hi the sum of the first i elements of b; 
T = (Ti , . . . , Tn): Ti the sum of the last n — 1 elements of b. 

Output. L: logical variable (L = False signals, that b is not graphical, while 
L = True shows, that the test could nor decide, whether b is graphical). 

Working variables, i: cycle variable; 
h: the actual value of hi; 

X = (Xi , X2, X3, X4, X5): Xj is the value of the actual Xij. 
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HEADSPLITTER-TEST(n, b, H, T, L) 

01 for i = 2 to n - 1 

02 h = [i/2J 

03 Xi = mm(HH, Tn — Ti, h(n — i)) 

04 X2 = mm(Hi - Hh, Tn - Tu [i - h) (n - i)) 

05 X3 = mm(h(t - h) 

06 X4 = Q) 



09 L = True 

10 return L 

The running time of the algorithm is 0(1) in best, and 0(n) in worst case. 

It is a substantial circumstance that the use of Lemma 3 and Lemma 4 
requires only linear time (while the earlier two theorems require quadratic 
time). But these improvements of Erdos-Gallai theorem decrease only the co- 
efficient of the quadratic member in the formula of the running time, the order 
of growth remains unchanged. 

Figure 2 contains the results of the running of Binomial-Test and Head- 
splitter-Test, further the values G(ti) and (the computation of the 

values of the function G(n) will be explained in Section 8). 

Figure 3 shows the relative frequency of the zerofree regular, binomial, head- 
splitted and graphical sequences compared to the number of regular sequences. 

3.4 Composite test 

Composite-Test uses approximate algorithms in the following order: Parity- 
Test, Binomial-Test, Positive-Test, Headsplitter-Test. 

COMPOSITE-TEST(n, b, L) 

01 PARITY-TEST(n, b,L) 

02 if L == False 

03 return L 

04 BiNOMiAL-TEST(n,b,H, L) 

05 if L == False 

06 return L 

07 HEADSPLITTER-TEST(n, b, H, T, L) 



07 
06 
07 
08 




if Hi > Xi + X2 + X3 + 2X4 + 2X5 or Xi + X2 > Ti 
L = False 
return L 
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08 if L == False 

09 return L 

10 L = True 

11 return L 

The running time of this composite algorithm is in all cases 0[ti). 

4 Properties of the approximate testing algorithms 

We investigate the efficiency of the approximate algorithms testing the regular 
algorithms. Figure 1 contains the number R(rL) of regular, the number E(rL) of 
even, and the number G(n) of graphical sequences for n = 1, . . . , 38. 

The relative efficiency of arbitrary testing algorithm A for sequences of given 
length n we define with the ratio of the number of accepted by A sequences of 
length n and the number of graphical sequences Gfn). This ratio as a function 
of TL will be noted by Xa(ti) and called the error function of A [34]. 

We investigate the following approximate algorithms, which are the compo- 
nents of Composite-Test: 

1) Parity-Test; 

2) Binomial-Test; 

3) Headsplitter-Test. 

According to (23) there are R(2) = 3 2-regular sequences: (1,1), (1,0) and 
(0,0). According to (25) among these sequences there are E(2) = 2 even se- 
quences. Binomial-Test accepts both even ones, therefore B(2) = 2. Both 
sequences are 2-graphical, therefore G(2) =2 and so the efficiency of Parity- 
Test (PT) and Binomial-Test (BT) is Xpt(2) = Xbt(2) =2/2 = 1, in this 
case both algorithms are optimal. 

The number of 3-regular sequences is R(3) = 10. From these sequences 
(2,2,2), (2,2,0), (2,1,1), (2,0,0) (1,1,0) and (0,0,0) are even, so E(3) = 
6. Binomial-Test excludes the sequences (2,2,0) and (2,0,0), so remains 
B(3) = 4. Since these sequences are 3-graphical, G(3) =4 implies Xpt(3) = | 
and Xbt(3) = 1. 

The number of 4-regular sequences equals to R(4) = 35. From these se- 
quences 16 is even, and the following 11 are 4-graphical: (3,3,3,3), (3,3,2,2), 
(3,2,2,1), (3,1,1,1,), (2,2,2,2), (2,2,2,0), (2,2,1,1), (2,1,1,0), (1,1,1,1), 
(1,1,0,0) and (0,0,0,0). From the 16 even sequences Binomial-Test also 
excludes the 5 sequences, so B(4) = G(4) = 11 and Xbt(4) = 1. 

According to these data in the case of n < 4 Binomial-Test recognizes all 
nongraphical sequences. Figure 2 shows, that for n < 5 we have B(n) = G(n), 
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Figure 2: Number of zerofree binomial, zerofree headsplitted and graphical 
sequences, further the ratio of the numbers of graphical sequences for neigh- 
bouring values of n 

that is Binomial-Test accepts the same number of sequences as the precise 
algorithms. If n > 5, then the error function of Binomial-Test is increasing: 
while Xbt(6) = (BT accepts one nongraphical sequence), Xbt(7) = 
(BT accepts 7 nongraphical sequences) etc. 
Figure 4 presents the average running time of the testing algorithms BT 
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Figure 3: The number of zerofree even sequences, further the ratio of the num- 
bers binomial/regular, headsphtted/regular and graphical/regular sequences 



and HT in secundum and in number of operations. The data contain the time 
and operations necessary for the generation of the sequences too. 
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Figure 4: Running time of Binomial-Test (BT) and Headsplitter-Test 
(HT) in secundum and as the number of operations for n = 1 , . . . , 20 

5 New precise algorithms 

In this section the zerofree algorithms, the shifting Havel-Hakimi, the parity 
checking Havel-Hakimi, the shortened Erdos-Gallai, the jumping Erdos-Gallai, 
the linear Erdos-Gallai and the quick Erdos-Gallai algorithms are presented. 

5.1 Zerofree algorithms 

Since the zeros at the and of the input sequences correspond to isolated ver- 
tices, so they have no influence on the quality of the sequence. This observation 
is exploited in the following assertion, in which p means the number of the 
positive elements of the input sequence. 

Corollary 5 If n > 1, the (bi,...,bn) n-regular sequence is n-graphical if 
and only i/ (bi , . . . , bp) is 'p- graphical. 

Proof. If all elements of b are positive (that is p = n), then the assertion 
is equivalent with Erdos-Gallai theorem. If b contains zero element (that is 
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p < n), then the assertion is the consequence of Havel-Hakimi and Erdos- 
Gallai algorithms, since the zero elements do not help in the pairing of the 
positive elements, but from the other side they have no own requirement. □ 

The algorithms based on this corollary are called Havel-Hakimi-Zerofree 
(HHZ), resp. Erdos-Gallai-Zerofree (EGZ). 

5.2 Shifting Havel-Hakimi algorithm 

The natural algorithmic equivalent of the original Havel-Hakimi theorem is 
called Havel-Hakimi Sorting (HHSo), since it requires the sorting of the 
reduced input sequence in every round. 

But it is possible to design such implementation, in which the reduction of 
the degrees is executed saving the monotonity of the sequence. Then we get 
Havel-Hakimi- Shifting (HHSh) algorithm. 

For the pseudocode of this algorithms see [37] . 

5.3 Parity checking Havel-Hakimi algorithm 

It is an interesting idea the join the application of the conditions of Erdos- 
Gallai and Havel-Hakimi theorems in such a manner, that we start with the 
parity checking of the input sequence, and only then use the recursive Havel- 
Hakimi method. 

For the pseudocode of the algorithm Havel-Hakimi-Parity (HHP) see 
[37]. 

5.4 Shortened Erdos-Gallai algorithm (EGSh) 

In the case of a regular sequence the maximal value of Hi is nfn— 1 ), therefore 
the inequality (2) certainly holds for i = n, therefore it is unnecessary to check. 

Even more useful observation is contained in the following assertion due to 
Tripathi and Vijai. 

Lemma 6 (Tripathi, Vijay [82]) If n > 1, then an n-regular sequence b = 
(bi . . . , bn) is n-graphical if and only if 

Hn even (13) 

and 

n 

Hi<min(Hi,i(i-1))+ ^ min(i,bfc) (i = 1 , 2, . . . , r), (14) 

k=i+l 
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where 

r = max (s I s(s — 1) < Hs) (15) 

1 <s<n 

Proof. If 1(1 — 1 ) > Hi, then the left side of (2) is nonpositive, therefore the 
inequahty holds, so the checking of the inequality is nonnecessary. □ 

The algorithm based on this assertion is called Erdos-Gallai-Shortened. 
For example if the input sequence is b = (5^"''), then Erdos-Gallai computes 
the right side of (2) 99 times, while Erdos-Gallai-Shortened only 6 times. 

5.5 Jumping Erdos-Gallai algorithm 

Contracting the repeated elements a regular sequence (bi , . . . , bn) can be writ- 
ten in the form (b^p...,b^), where bi, < ••• < b^^, ei, Cq > 1 and 

ei H h eq = n. Let = ei H \- Cj (j = 1, . . . , q). 

The element bt is called the checking points of the sequence b, if i = n or 
1 <i<n— 1 esbi> b^+i . Then the checking points are bg, , . . . , bg^ . 

Theorem 7 (Tripathi, Vijay [82]) An n-regular sequence b = (bi, . . . ,bn) is 
n-graphical if and only if 

Hn even (16) 

and 

n 

Hg. -gi(gi-l) < Y_ min(i,bk) (i = l,...,q). (17) 

k=gi+l 

Proof. See [82]. □ 

Later in algorithm Erdos-Gallai-Enumerating we wiU exploit, that in 
the inequality (17) gq is always n, therefore it is enough to check the inequality 
only up to 1 = q — 1 . 

The next program implements a quick version of Erdos-Gallai algorithm, 
exploiting Corollary 5, Lemma 6 and Lemma 7. In this paper we use the 
pseudocode style proposed in [17]. 

Input, tl: number of vertices (n > 1); 
b = (bi , . . . , bn): an n-even sequence. 

Output. L: logical variable (L = False signalizes, that, b is not graphical, 
while L = True shows, that b is graphical). 

Working variables, i and j: cycle variables; 
H = (Ho, Hi , . . . , Hn): H^ is the sum of the first i elements of b; 
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C: the degree capacity of the actual tail; 

bn+i : auxiliary variable helping to decide, whether bn is a jumping element. 
Erdos-Gallai- JuMPiNG(n, b, H, L) 

01 Hi = bi // hues 01-06: test of parity 

02 for t = 2 to n 

03 Hi = Hi_i + bi 

04 if Hn odd 

05 L = False 

06 return L 

07 bn+i = — 1 // lines 07-20: test of the request of the head 

08 i = 1 

09 while t < n and i(i — 1 ) < Hi 

10 while bi == bi+i 

11 i = i+l 

12 C = 

13 for j = i + 1 to n 

14 C = C + min(j,bj) 

15 if Hi > t(t- 1) + C 

16 L = False 

17 return L 

18 1 = 1+1 

19 L = True 

20 return L 

The running time of EGJ varies between the best 0(1 ) and the worst &[n^]. 
5.6 Linear Erdos-Gallai algorithm 

Recently we could improve Erdos-Gallai algorithm [35, 37]. The new al- 
gorithm Erdos-Gallai-Linear exploits, that b is monotone. It determines 
the capacities Ci in constant time. The base of the quick computation is the 
sequence w(b) containing the weight points wi of the elements of the input 
sequence b. 

For given sequence b let w(b) = (wi, . . . ,Wn-^ ), where Wi gives the index 
of bi^ having the maximal index among such elements of b which are greater 
or equal to i. 

Theorem 8 (Ivanyi, Lucz [35], Ivanyi, Lucz, Mori, Soter [37]) Ifn> 1, then 
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the n-regular sequence (bi , . . . , bn) is n-graphical if and only if 

Hn is even (18) 

and if i > Wi, then 

Hi<l(l-1)+Hn-Ht, 

further if i < Wi, then 

Hi < i(i - 1 ) + i[wi - I) + Hn - H^^. 

Proof. (18) is the same as (1). 

During the testing of the elements of b by Erdos-Gallai-Linear there 
are two cases: 

• if i > Wi, then the contribution Ct = ^{J^^t+i min(i, b^) of the tail of b 
equals to Hn — H^, since the contribution Ck of the element is only b]^. 

• if i < wi , then the contribution of the tail of b consists of contributions of 
two types: Ci+i , . . . , c^^ are equal to i, while Cj = bj for j = + 1 , . . . , n. 

Therefore in the case n — 1 > t > Wi we have 

Q = Hn - Hi, (19) 

and in the case 1 < i < Wi 

Ci = i(wi-i) + Hn-H^^. (20) 

□ 

The following program is based on Theorem 8. It decides on arbitrary n- 
regular sequence whether it is n-graphical or not. 

Input, n: number of vertices (n > 1); 
b = (bi , . . . , bn): n-regular sequence. 

Output. L: logical variable, whose value is True, if the input is graphical, 
and it is FALSE, if the input is not graphical. 

Work variables, i and j: cycle variables; 
H = (Hi , . . . , Hn): Hi is the sum of the first i elements of the tested b; 
bo: auxiliary element of the vector b 

w = (wi , . . . , Wn-i ): Wi is the weight point of bi, that is the maximum of the 
indices of such elements of b, which are not smaller than i; 
Ho = 0: help variable to compute the other elements of the sequence H; 
bo = n — 1 : help variable to compute the elements of the sequence w. 
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ERD6s-GALLAI-LlNEAR(rL, b, L) 

01 Ho = // line 01: initialization 

02 for t = 1 to n // lines 02-03: computation of the elements of H 

03 Hi = Hi_i + bi 

04 if Hn odd // lines 04-06: test of the parity 

05 L = False 

06 return L 

07 bo = TL — 1 // line 07: initialization of a working variable 

08 for t = 1 to n // lines 08-12: computation of the weights 

09 if bi < bi_i 

10 for j = bi 1 downto bi + 1 

11 Wj = i — 1 

12 w^. = i 

13 for j = bn downto 1 / / lines 13-14: large weights 

14 Wj = n 

15 for t = 1 to n // lines 15-23: test of the elements of b 

16 if I < Wi // lines 16-19: test of indices for large Wi's 

17 if Hi > i(i - 1 ) + l(wi - I) + Hn - H^. 

18 L = False 

19 return L 

20 if I > Wi // lines 20-23: test of indices for small Wi's 

21 if Hi >i(i-1) +Hn-Hi 

22 L = False 

23 return L 

24 L = True / / lines 24-25: the program ends with the value True 

25 return L 



Theorem 9 (Ivanyi, Lucz [35], Ivanyi, Lucz, Mori, Soter [37]) Algorithm 
Erdos-Gallai-Linear decides in Q[n) time, whether an n-regular sequence 
b = (bi , . . . , bn) is graphical or not. 

Proof. Line 1 requires 0(1) time, lines 2-3 @{r\.) time, lines 4-6 0(1) time, 
line 07 0(1) time, lines 08-12 0(1) time, lines 13-14 0(n) time, lines 15-23 
0(n) time and lines 24-25 0(1) time, therefore the total time requirement of 
the algorithm is 0(n). □ 

Since in the case of a graphical sequence all elements of the investigated 
sequence are to be tested, in the case of RAM model of computations [17] 
Erdos-Gallai-Linear is asymptotically optimal. 
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6 Running time of the precise testing algorithms 

We tested the precise algorithms determining their total running time for all 
the even sequences. The set of the even sequences is the smallest such set of 
sequences, whose the cardinality we know exact and explicite formula. The 
number of n-bounded sequences K(rL) is also known, but this function grows 
too quickly when n grows. 

If we would know the average running time of the bounded sequences we 
would take into account that is is sufficient to weight the running times of the 
regular sequences with the corresponding frequencies. For example a homoge- 
neous sequence consisting of identical elements would get a unit weight since 
it corresponds to only one bounded sequence, while a rainbow sequence con- 
sisting is n different elements as e.g. the sequence n, n — 1 , . . . , 1 corresponds 
to n! different bounded sequences and therefore would get a corresponding 
weight equal to n!. 

We follow two ways of the decreasing of the running time of the precise 
algorithms. The first way is the decreasing of the number of the executable 
operations. The second way is, that we try to use quick (linear time) prepro- 
cessing algorithms for the filtering of the sequences in order to decrease of the 
part of sequences requiring the relative slow precise algorithms. 

For the first type of decrease of the expected running time is the short- 
ening of the sequences and the application of the checking points, while for 
the the second type are examples the completion of HH algorithm with the 
parity checking or the completion of the EG algorithm with the binomial and 
headsplitted algorithms. 

In this section we investigate the following precise algorithms: 

1) Havel-Hakimi-Shorting (HHSo). 

2) Havel-Hakimi-Shifting (HHSh). 

3) Erdos-Gallai algorithm (EG). 

4) Erdos-Gallai-Jumping algorithm (EGJ). 

5) Erdos-Gallai-Linear algorithm (EGL). 

Figure 5 contains the total number of operations of the algorithms HHSo, 
HHSh, EG, and EGL required for the testing of all even sequences of length 
n = 1 , . . . , 15. The operations necessary to generate the sequences are included. 

Gomparison of the first two columns shows that algorithm HHSh is much 
quicker than HHSo, especially if n increases. Gomparison of the third and 
fourth columns shows that we get substantial decrease of the running time 
if we have to test the input sequence only in the check points. Finally the 
comparison of the third and fifth columns demonstrates the advantages of a 
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Figure 5: Total number of operations as the function of n for precise algorithms 
HHSo, HHSh, EG, EGJ, and EGL. 
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Figure 6: Total and amortized running time of Erdos-Gallai-Linear in 
secundum, resp. in the number of executed operations 
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Figure 7: Distribution of the even nongraphical sequences according to the 
number of tests made by Erdos-Gallai-Jumping to exclude the given se- 
quence for n = 3, 15 

hnear algorithm over a quadratic one. 

Figure 6 shows the running time of Erdos-Gallai-Linear in secundum 
and operation, and also the amortized number of operation/even sequence. 

The most interesting data of Figure 6 are in the last column: they show 
that the number of operations/investigated sequence/length of the investigated 
sequence is monotone decreasing (see [69]). 

Figure 7 shows the distribution of the E(tl) — G(n) even nongraphical se- 
quences according to the number of tests made by Erdos-Gallai-Jumping 
to exclude the given sequence for n = 3, . . . , 15 vertices, idn) = fi gives the 
frequency of even nongraphical sequences of length n, which requeired exactly 
i round of the test. 

These data show, that the maximal number of tests is about j in all lines. 

Figure 8 shows the average number of required rounds for the nongraphical, 
graphical and all even sequences. The data of the column belonging to G(rL) 
are computed using Lemma 17. It is remarkable that the sequences of the 
coefficients are monotone decreasing in the last three columns. 

Figure 9 presents the distribution of the graphical sequences according to 
their first element. These data help at the design of the algorithm Erdos- 
Gallai-Enumerating which computes the new values of G (n) (in the slicing 
of the computations belonging to a given value of n). 
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Figure 8: Weighted average number of tests made by Erdos-Gallai- Jumping 
while investigating the even sequences for n = 3, . . . , 15 
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Figure 9: The distribution of the graphical sequences according to bi for n = 
1, 12 



We see in Figure 9 that from n = 6 the multiplicities increase up to n — 2, 
and the last positive value is smaller then the last but one element. 
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7 Enumerative results 

Until now for example Avis and Fukuda [4], Barnes and Savage [5, 6], Burns 
[13], Erdos and Moser [59], Erdos and Richmond [22], Frank, Savage and Sel- 
ers [25], Kleitman and Winston [47], Kleitman and Wang [46], Metropolis and 
Stein [56], R0dseth et al. [68], Ruskey et al. [69], Stanley [78], Simion [71] and 
Winston and Kleitman [86] published results connected with the enumeration 
of degree sequences. Results connected with the number of sequences investi- 
gated by us can be found in the books of Sloane es Ploffe [76], further Stanley 
[79] and in the free online database On-line Encyclopedia of Integer Sequences 
[73, 74, 75] 

It is easy to show, that if I, u and ra are integers, further u > I, m > 1 , and 
I < < u for i = 1 , . . . , ra, then the number of (I, u, m)-bounded sequences 
a = (ai , . . . , Qra) of integer numbers K(l, u, m) is 

K(l,u,ra) = (u-l + 1)^. (21) 

It is known (e.g. see [39, page 65]), that if I, u and m are integers, further 
u > I and m > 1 , and u > bi > • • • > bn > 1, then the number of (I, u, m)- 
regular sequences of integer numbers R(l, u, m) is 

R(i,u,m)=("-;;^+'"). (22) 

The following two special cases of (22) are useful in the design of the algo- 
rithm Erdos-Gallai-Enumerating. 

If n > 1 is an integer, then the number of R(0,tl— l,n)-regular sequences is 

mn-^,n)=R[n) = (^'^~'^y (23) 

If n > 1 is an integer, then the number of R(l , n — 1 , n)-regular sequences is 

R(l,n-l,n] =RJn) = (^"^ ~ ^) • (24) 

In 1987 Ascher derived the following explicit formula for the number of 
n-even sequences E(n). 



Lemma 10 (Ascher [3], Sloane, Pfoffe [76]) // 
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Lemma 11 lemma-En n > 1, then the number of n- even sequences E(n) is 

Proof. See [3, 76]. □ 

At the designing and analysis of the resuhs of the simulation experiments 
is useful, if we know some features of the functions Rfnl and Efnl. 



Lemma 12 // n > 1 , then 



R(n + 2) R(n+1) 
R(n+1)^ R(n) ' 

(27) 



n— >oo 



Rfu) 



further 



4^ / 1 \ f 1 

1-^ <W<^ 1-^-^)- (28) 



yjAnn V 2ny V^tttl V 8n + 8 

Proof. On the base of (23) we have 

R(n + 2) ^ (2n + 3)!(n + 1)n! ^ 4n + 6 ^ ^_ 
R(n+1) (n + 2)!(n+l)!(2n + l)! n + 2 n + 2' ^ ^ 

from where we get directly (26) and (27). □ 
Using Lemma 13 we can give the precise asymptotic order of growth of E(u). 



Lemma 13 //n > 1, then 



further 



^(" + ^^ > ^^" + 1) (30) 



lim%±il=4, (31) 
n^oo En) 



'l-D3(n))<E(n)<^(l-D4(n)), (32) 



V7rn \/Tir\. 
where D3(n) and D4(n) are monotone decreasing functions tending to zero. 
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Proof. The proof is similar to the proof of Lemma 12. □ 

Comparison of (23) and Lemma 13 shows, that the order of growth of num- 
bers of even and odd sequences is the same, but there are more even se- 
quences than odd. Figure 1 contains the values of Rln), E(n) and E(n)/R(rL) 
for n = 1 , . . . , 37. 

As the next assertion and Figure 1 show, the sequence of the ratios E(u)/R(n) 
is monotone decreasing and tends to j. 



Corollary 14 // n > 1 , then 



< 



Efnl 



and 



H(n + i: 

R(n-h 1) ^ R(rL) 
Efnl 1 



lim 
n->oo Kin] 



(33) 
(34) 



Proof. This assertion is a direct consequence of (23) and (25). 



□ 



The expected value of the number of jumping elements has a substantial in- 
fluence on the running time of algorithms using the jumping elements. There- 
fore the following two assertions are useful. 

The number of different elements in an n-bounded sequence b is called the 
rainbow number of the sequence, and it will be denoted by Tn(b). 

Lemma 15 Letb be a random n-bounded sequence. Then the expectation and 
variance of its rainbow number are as follows. 



E[rn(b)] =n 



1 



1 



n 



Var[rn(b)] = n 1 



n 



1 



1 



1 



+ 0(1), 



(35) 



n 



-Fnfn- V 



e V e. 



n 



2n 



(36) 



Proof. Let E,i denote the indicator of the event that number i is not contained 
in a random n-bounded sequence. Then the rainbow number of a random 
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sequence is 'n — Y_^=o hence its expectation equals ti — Y.u^o ^[^i]- Clearly, 



n 



holds independently of i, thus 





1 - 




n 











On the other hand, 
Var[rn(b)] = Var 

Here 



n-l ■ 
.i=0 . 



n-1 



}^Var[y+2 Y. cov[^t,£,j: 

0<i<j<n-l 



Var[^t] 



i=0 



n 



and 



cov[£,i,L] =E[£,iL]-E[£,i]E[L] = 1 



n 

2\ n 

n, 



n 



2n 



(37) 
(38) 

(39) 

(40) 
(41) 



implying (36). □ 

We remark that this lemma answers a question of Imre Katai [40] posed in 
connection with the speed of computers having interleaved memory and with 
checking algorithms of some puzzles (e.g sudoku). 

Lemma 16 The number of (0,n — ] ,m)-regular sequences composed from k 
distinct numbers is 



m- 1 



, k = 1 , . . . , n. 



(42) 



In other words, the distribution of the rainbow number r-r^{h] of a random 
(0, n — 1 , rrCj-regular sequence b is hypergeometric with parameters n + m — 1 , 
n, and m. 

Proof. The k-set of distinct elements of the sequence can be selected from 
{0, 1 , . . . , n — 1} in (]J) ways. Having this values selected we can tell their 
multiplicities in {^Z-l) ways. Let us consider the k blocks of identical elements. 
The first one starts with bi , and the starting position of the other k — 1 blocks 
can be selected in {'^_^) ways. □ 

From this the expectation and the variance of a random n-regular sequence 
follow immediately. 



On Erdos-Gallai and Havel-Hakiini algorithms 



257 



Corollary 17 Let b be a random n-regular sequence. Then the expectation 
and the variance of its rainbow number rn[b] are as follows: 

TL^ TL TL TL 

= fe? = i + I28n^ - I28n + 32 = ¥ + 

Lemma 18 Lei h be a random n-regular sequence. Let us write it in the form 
b = (b^^ , . . . ,b^''). Then the expected value of the exponents Cj is 

E[ej |r(b) > j] =4 + 0(1). (45) 

Proof. Let c(n, j) denote the number of n-regular sequences with rainbow 
number not less than j. By Lemma 16, 



Let us turn to the number of n-regular sequences with rainbow number not 
less than j and ej = £. This is equal to the number of (0, n— 1 , n—l + ^ )-regular 
sequences containing at least j different numbers, that is, 

£0 (:::)• 

Prom this the sum of ej over all n-regular sequences with ej > is equal to 

{=1 k=j ^ ^ ^ ^ k=j ^ ^ f=i ^ ^ ^ 

n 



k=j 

This can also be seen in a more direct way. Consider an arbitrary n-regular 
sequence with at least j + 1 blocks, then substitute the elements of the j + 1st 
block with the number in the jth block (that is, concatenate this two adjacent 
blocks) and delete one element from the united block; finally, decrease by 1 
all elements in the subsequent blocks. In this way one obtains an n-regular 
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sequence with at least j blocks, and it easy to see that every such sequence is 
obtained exactly ej times. 

Thus the expectation to be computed is just 



C(T1+I,j + 1) 

c(rL,j) 

Clearly, c(n, 1 ) = R(0, n — 1 , n) = ^ j , hence 



2ti — 1 \ t-^ /n\ /n — 1 \ /2tl — 1 



(49) 



as TL — > oo. This is asymptotically equal to 



2rL + r 
n + 1 / 4n + 2 



4 ^=4 + 0(1). (51) 



2rL-l\ TL+l n+1 
n J 

□ 

It is interesting to observe that by (43) the average block length in a random 
n-regular sequence is 

-Z^j = 4T~2 (52) 
r ■pj- r[bj 

approximately, as n — ) oo. This fact could be interpreted as if blocks in the 
beginning of the sequence were significantly longer. However, fixing rn(b) = r 
we find that the lengths of the r blocks are exchangeable random variables with 
equal expectation n/r. At first sight this two facts seem to be in contradiction. 
The explanation is that exchangability only holds conditionally. Blocks in the 
beginning do exist even for smaller rainbow numbers, when the average block 
length is big, while blocks with large index can only appear when there are 
many short blocks in the sequence. 

The following assertion gives the number of zerofree sequences and the ratio 
of the numbers of zerofree and regular sequences. 

Corollary 19 The number of the zerofree n-regular sequences Rz;(n) is 

Rz(n) = (^^~^) (53) 
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and 

li- ¥^ = (54) 

Proof. (53) identical with (22), (54) is a direct consequence of (22) and (23). 

□ 

As the experimental data in Figure 3 show, ~ \- 
The following lemma allows that the algorithm Erdos-Gallai-Enumera- 
ING tests only the zerofree even sequences instead of the even sequences. 

Lemma 20 7/ n > 2, then the number of the n-graphical sequences G (n) is 

G(n) = G(n-l) + GJn). (55) 

Proof. If an n-graphical sequence b contains at least one zero, that is = 0, 
then b' = (bi , . . . , b^-i ) is (n — 1 )-graphical or not. If a = (ai , . . . , Un-i ) is 
an (n — 1 )-graphical sequence, then a' = (ai , . . . , an-i , 0) is n-graphical. 

The set of the n-graphical sequences S consists of two subsets: set of zerofree 
sequences Si and the set of sequences S2 containing at least one zero. There 
is a bijection between the set of the (n — 1 )-graphical sequences and such n- 
graphical sequences, which contain at least one zero. Therefore |S| = |Si I+IS2I = 
GJn) + G(n-l). □ 

Corollary 21 //n > 1, then 

n 

G(n) = 1 +_^GJn). (56) 

Proof. Concrete calculation gives G(l) = 1. Then using (55) and induction 
we get (56). □ 
A promising direction of researches connected with the characterization of 
the function G(n) is the decomposition of the even integers into members 
and the investigation, which decompositions represent a graphical sequence 
[5, 6, 13]. Using this approach Burns proved the following asymptotic bounds 
in 2007. 

Theorem 22 (Burns [13]) There exist such positive constants c and C, that 
the following bounds of the function G(n) is true: 

— < G(n) < (57) 

cn [lognj'-V'"- 
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Proof. See [13]. □ 
This result implies that the asymptotic density of the graphical sequences 
is zero among the even sequences. 

Corollary 23 // n > 1 , then there exists a positive constant C such that 

eR < (bi^ ^^^^ 

and 

lim ^=0. (59) 

Proof. (58) is a direct consequence of (25) and (58), and (58) implies (59). □ 
As Figure 1 and Figure 3 show, the convergence of the ratio G(u)/E(n) is 
relative slow. 



8 Number of graphical sequences 

Erdos-Gallai-Enumerating algorithm (EGE) [37] generates and tests for 
given n every zerofree even sequence. Its input is n and output is the number 
of corresponding zerofree graphical sequences Gz(rL). 

The algorithm is based on Erdos-Gallai-Linear algorithm. It generates 
and tests only the zerofree even sequences, that is according to Corollary 5 
and Figure 3 about the 25 percent of the n-regular sequences. 

EGE tests the input sequences only in the checking points. Corollary 17 
shows that about the half of the elements of the input sequences are check 
points. 

Figure 3 contains data showing that EGE investigates even less than the 
half of the elements of the input sequences. 

Important property of EGE is that it solves in 0(1) expected time 

1. the generation of one input sequence; 

2. the updating of the vector H; 

3. the updating of the vector of checking points; 

4. the updating of the vector of the weight points. 

We implemented the parallel version of EGE (EGEP). It was run on about 
200 PC's containing about 700 cores. The total running time of EGEP is 
contained in Figure 10 

The pseudocode of the algorithm see in [37]. The amortized running time of 
this algorithm for a sequence is 0(1), so the total running time of the whole 
program is 0(E(n)). 
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TL 


running time (in days) 


number of slices 


24 


7 


415 


25 


26 


415 


26 


70 


435 


27 


316 


435 


28 


1130 


2001 


29 


6733 


15119 



Figure 10: The runnng time of EGEP for n = 24, . . . , 29 

9 Summary 

In Figure 1 the values of R(n.) up to n = 24 are the elements of the sequence 
A001700 of OEIS [73], the values of E(tl) up to n = 23 are the elements of the 
sequence A005654 [75] of the OEIS, and in Figure 2 the values G(rL) are up to 
n = 23 are the elements of sequence A0004251-es [74] of OEIS. The remaining 
values are new [37, 36]. 

Figure 2 contains the number of graphical sequences G (n) for n = 1 , . . . , 29, 
and also G (n + 1 )/G (n) for n = 1 , . . . , 28. 

The referenced manuscripts, programs and further simulation results can be 
found at the homepage of the authors, among others at http: //compalg. inf . elte.hu/~tony/Kutatas/Ei 
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